SIGNIFICANCE  AND  EXPLANATION 


Serial  systems  arise  naturally  in  practice  in  engineering 
and  physics.  Therefore  it  is  of  substantial  significance  to  be 
able  to  efficiently  utilize  data  obtained  on  individual  com¬ 
ponents  for  the  purpose  of  obtaining  an  overall  assessment  of 
the  reliability  of  the  system.  The  methods  of  this  paper  may 
be  employed  for  this  purpose. 
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.  1.  Introduction  and  Summary 

A  problem  of  f undamental  interest  to  practitioners  in 
reliability  is  the  statistical  estimation  of  the  reliability  of  a 
system  using  experimental  data  collected  on  subsystems.  In  this 
paper,  the  subsystem  data  available  consists  of  a  sequence  of 
Bernoulli  trials  in  which  a  "one"  is  recorded  if  the  subsystem 
functions  and  a  zero  is  recorded  if  the  subsystem  falls.  Thus 
for  each  of  the  k  subsystems  composing  the  system,  the  data  pro¬ 
vided  consists  of  the  pair  (n^Y^),  i«l,2,...k,  where  Y±  is 
binomially  distributed  (n^p^.  We  assume  that  Y^,Y2,...,Yfc  are 
mutually  independent  random  variables. 

The  magnitude  of  interest  in  this  problem  is  easily  evidenced 
by  the  extensive  literature  devoted  to  it.  In  this  regard,  see 
the  survey  paper  by  Harris  (1977)  and  Section  10.4  of  the  book  by 
Mann,  Schafer,  and  Singpurwalla  (1974).  In  addition,  the  Defense 
Advanced  Research  Projects  Agency  has  recently  Issued  a  Handbook 
for  the  Calculation  of  Lower  Statistical  Confidence  Bounds  on 
System  Reliability  (1980). 

Historically,  the  first  significant  work  on  this  problem  was 
produced  by  Buehler  (1957).  However,  Buehler's  method  as  des¬ 
cribed  in  that  paper  is  difficult  to  implement  computationally 
when  k>2. 
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We  proceed  by  describing  Buehler's  method  in  Section  2.  In 
Section  3  we  specialize  to  series  systeas,  that  is,  a  system  which 


fails  whenever  at  least  one  subsystem  fails.  Sudakov's  (1974) 
results  are  extended  in  Section  4  and  employed  to  exhibit  some 
optimality  properties  of  the  L ind s t r om-Madd en  method  (see  Lloyd 
and  Lipow  (1962))  for  constructing  lower  confidence  bounds  for 
the  reliability  of  series  systems  of  stochastically  independent 
subsystems.  Some  numerical  examples  are  given  in  Section  5  and 
the  results  needed  for  this  generalization  of  Sudakov's  Theorem 
are  provided  in  the  Appendix  to  this  paper. 

2.  Buehler's  Method  for  Lower  Confidence  Bounds 

A  system  composed  of  k  independent  subsystems  is  said  to  be 
a  coherent  system  (with  respect  to  the  specified  decomposition 
into  subsystems),  if  the  system  fails  when  all  subsystems  fail 
and  the  system  functions  when  all  subsystems  function;  and  replac¬ 
ing  a  defective  subsystem  by  a  functioning  subsystem  can  not 
cause  a  functioning  system  to  fail.  Coherent  systems  are  des¬ 
cribed  in  Birnbaum,  Esary  and  Saunders  (1961)  and  Barlow  and 
Proschan  (1975). 

To  any  system  one  can  associate  a  function,  h(p)  = 

h(pl » p2 ’ *  * ’ ’ Pp) *  iel»2,...»k,  where  h(p)  is  the  reliabil- 

t  h 

ity  of  the  system  when  pi  is  the  probability  that  the  i— 
subsystem  functions.  It  is  well-known  that  if  the  system  is 
coherent , 

0  <  h(p)  <  1  , 

h(0 . 0)  -  0,  h (1 ,  .  .  .  ,  1 )  «=  1  , 
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7  . 


and  h(p^,...,p^)  is  non-decreasing  in  each  variable 


For  coherent  systems,  Buehler's  method  may  be  described  as 

follows:  The  observed  outcome  (y.,.,.,y.)  can  assume  any  of 

k  lk 

N  ■  II  (n.+l)  values,  since  y.  »  0,1,  ...,n..  For  convenience, 
1-1  1  1  1 
we  denote  n^-y^  by  x^»  1*1,2, ...,k. 

A  partition  '(Aj'tAg,  •••  ,A  )  ,  s>l,  of  the  N  possible  outcomes 
Is  said  to  be  a  monotonic  partition,  that  is,  A1<A2<...<Ag  if 

*°’° . e  Al*  (ni*n2»  *  •  *  *nk^  £  As  and  if  xi  "  ^xn . xlk*’ 

x,  -  (x,,  ,  . . .  ,xolf )  with  x^  <  x^,  i*l,2 . k,  then  x,  C  A4 


implies  C  A ^ ,  j  >  1. 


k  fn.1}  n  .-x  x.  k  n 
f(x;5)  *  P-(X*x)  -  n  p  1  \  1  -  It 

p  1-lUJ  i-ily 


*  rn  i  rt  -j-yj 
h  I  IP±  q± 


(2.1) 


and  for  1  <  n  <  s-1,  let 


a  *  inf 
n 


jh(p) I  l  f (x. ;p)  -  a) 
1  |xieAi,l<n  1  > 


(2.2) 


and  a  -  0. 
s 

Each  such  partition  may  be  identified  with  a  function  defined  on 
the  set  of  sample  outcomes  by  defining  the  ordering  function 
g(x),  where 


g(x)  -  n  if  x  e  A  ,  l<n<s  ; 

n  —  — 


(2.3) 


!■  1 


obviously  g(x)  inherits  the  monotonicity  properties  of  the 
partition. 

Subsequently  it  will  be  convenient  to  use  ordering  functions 
g ( x )  such  that  the  range  of  g(x)  will  be  a  finite  set  of  real 
numbers,  r ^<r . . . <r g .  With  no  loss  of  generality,  we  can  identify 


the  se 


ts  A±  by  defining  A^  -  jx|g(x)  -  r^,  i-1,2, ...,s. 


We  can 


now  establish  the  following  theorems. 
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Theorem  2.1.  Let  x  be  distributed  by  (2.1).  Then  a  „  Is  a 
-  g(x> 

(l-o)  lower  confidence  bound  for  h(p).  If  b  is  also  a  (1-a) 

g  (X  ) 

lower  confidence  bound  for  h(p),  then  b^<a^,  l<i<s. 

Proof :  Fix  p  and  let  n(p)  be  the  smallest  integer  such  that 


(  n (p)  ) 

?\X  £  ±Vml  Ai)  -  °  * 


(2.4) 


>-{x  e  U  A  }  >  1-a  . 
pl  i-n(p)  ij 


(2.5) 


Dn  -  «  JA,}  »  a} 


(2.6) 


Then  D  is  a  1-a  confidence  set  for  p,  since 

Pp{?  e  °g(X)}  "  Pp{g(*>  -  n(p)}  -  1-01  *  (2*7) 

This  establishes  the  first  part  of  the  conclusion.  Further,  since 

h(p)  is  continuous  and  0<p^<l,  the  infimum  in  (2.2)  is  attained. 

Now  assume  that  1.  is  the  smallest  index  such  that  b  >a  , 

11 

l<ij<8-l.  Then,  for  some  pQ,  p^, 

b  >  inf{h(p)  |  l  f(x;p)  *=  a)  -  h(p  )  » 

1  1  x1£A1,i<ii  J  ° 


I  f (x;p  )  >  a,  h(p. )  <  b, 

1  1 


xieA1,i<i1 


Therefore 


<  bxC*>}  i  5  >  - 


a  contradiction. 


■fl-h(p)  l  f  (x.  ;p)  -  a}. 

£  xicAi*1-n  ' 

(1-a)  upper  confidence  bound  for  l-h(p),  the  unreliability 


Remark.  Let  d  «  sup 
-  n 


Then  d  is  a 
n 


Let  A 


{ 


xeE^,  1“1 1 


,2,...,kj  and  let  g(x)  be  contin¬ 
uous  on  A  (the  closure  of  A)  and  strictly  increasing  in  each 
variable  for  xeA.  g(x)  is  to  be  regarded  as  an  ordering  function 
as  described  immediately  preceding  Theorem  2.1.  We  require  the 
following  additional  property  of  g(x). 

Fix  x  £A.  Let  g(x  )  <  g (a . , 0 , . , . , 0)  =  g. .  Then 

O  OX  X 

g (y J , 0, . . . , 0)  =  g(xQ)  has  a  unique  solution  in  y^  .  Proceeding 
recursively,  let  i^  <  y^  and  define  y^  ■  y2^1^  aS  t*'e  so^-ut*on 
of  g(xQ)  -  g(i1,y2,0,...,0).  For  each  l<jfk  and  ij_i  5  yj-l’ 


lj-2  -  yj  -  2  ’  *  *  *  *  il  -  yl* 


let  y 


i 


yj  ^tl»i2 . 1j~l)  the 


solution  of 


g(xo)  “  g (i^ * ^2  *  *  *  *  * ^ J -1  * y j * ^ *  *  *  *  * ^  *  (2.8) 

We  require  that  the  equations  indicated  in  (2.8)  have 
unique  solutions  for  each  yj . 

Then  define 

£ y^]  [y2J  tykl 

F (x  ;p)  “J  l  •••  Z  f(i;p)  ,  (2.9) 

V0  v°  v° 

where,  for  j>l,  y^  *  y^  (i^,  i2  , . .  , ,  i^  .  Let 

f*(x  ;a)  -  sup  F(x  ;p)  ,  0<a<l  .  (2.10) 

°  h(p)-a 

Then  we  have 

Theorem  2.2.  If  x  satisfies  inf  f*(x  ;a)  **  0,  sup  f*(x  ;a)  *  1 

°  0<a<l  0<a< 1  0 

it 

and  f  (xQ;a)  is  a  strictly  increasing  function  of  a,  and  if 

x  e  A„  where  g(x)  determines  (A. , A- , . . . , A  ) ,  and  if 
on  or  x  i.  s 
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b  ■  inf 


(h(p) I  l  f  (x  ,?)  =  a)  , 
|xeA.,i<n  1  J 


(2.11) 


t  hen  we  have 


f*(xQ;b)  =  a 


Proof:  Since  the  infimun  in  (2.11)  is  attained,  there  is  a  pQ 

*«»  m*  A  ^ 

such  that  b  =  h(p  )  and  F(x  ;p  )  s  a.  Then  £  (x  ,  b)  >  a.  If 
vro  oro  o  - 

f  (x^.b)  >  a,  there  exists  p^  ,  with  a  =  h(p^),  a<b  and 

*  -w 

f  (xQ;a)  =  a  contradicting  (2.11). 

Obviously,  the  above  discussion  can  easily  be  modified  to 
obtain  upper  confidence  bounds  on  the  unreliability  l-h(p)  by 
replacing  inf  by  sup  in  (2.11)  and  requiring  that  f  (xQ;a)  be  a 
strictly  decreasing  function  of  a,  0<a<l. 


3.  Applications  to  Series  Systems 

k 

For  a  series  system  h(p)  =  H  p..  Further,  throughout  this 

1=1  1 

section  we  assume  that  g(x)  satisfies  the  conditions  necessary  to 

insure  that  the  solutions  for  y^,...,y^  indicated  in  (2.8)  are 

unique.  Then  we  have  the  following  theorem. 

k 

Theorem  3.1.  If  h(p)  «  II  p  ,  then  inf  f  (x  ;a)  =  0, 

*  _  iKl  0<a<l  ° 

sup  f  (x  ,a)  **  1  and  f*(x  ;a)  is  strictly  increasing  in  a, 

0<a<l  ° 

whenever  xQ  *>  (xq1  ,  .  .  .  ,  satisfies  xqj  <  n^  ,  j«l,2,...,k. 

Proof .  Since  h(p)  *  1  if  and  only  if  p  ^  *  1,  i=l,2,...,k,  it 


follows  from  (2.1)  that 


lim  sup  F(x  ;p)  ■  1  . 
a-+l  h(p)**a  w 
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Similarly,  h(p)  *  0  if  and  only  if  at  least  one  p^  ■  0, 

i-1,2 . k.  Since  F(xo;p)  <  P-jx^n^  «  l-P-jx^n^  =  1-q^, 

we  have 


lim  sup  F  C x  J  p )  ■  0  . 
a-*0  h(p)Ka  c 


To  show  that  f  (xq ; a)  is  strictly  increasing  in  a,  consider 

A  ^  _ 

0<a<b<l  and  let  pa  -  (pal » *  *  * ‘ pak^  satisfy  f  (*o;a>  “  F(Xo;Fa^‘ 

•ft 

Similarly,  let  satisfy  f  (xo;b)  -  F(x0;pb^‘  Let 
1  *  |ij , ij » • • • » ir |  be  any  non-empty  set  of  indices  such  that 

Pa.  (~)1//r  <  1  and  let  IC  be  the  remaining  Indices.  Then 


(  n  p 

JCI 


)  n 

jclc 


b  . 


(3.1) 


From  the  monotone  likelihood  ratio  property  of  the  binomial 
distribution, 

F(xo;pa)  <  F(£o;p*)  , 

where  the  components  of  p*  are  given  by  (3.1).  Then 

F(xo;p*)  <  sup  F(xQ:p)  -  F(xo;rb)  “  f* (xQ;b)  . 
h  (p)*=b 


4.  Sudakov's  Methoc 


Let 

Vr,s)  ■  ¥77777  C  ' 

Then  if  y  is  an  integer,  y<n,  we  have 


][  (?)  Pn”1q1  -  I  (n-y,y+l)  . 

i-0  1  P 
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For  0<y<n,  real,  define  u(n,y,a)  by  a  -  1u(n,y ,o) ^n“y »y+1^ • 

Thus,  for  Integer  values  of  y,  u(n,y»ct)  is  a  lOO(l-a)  percent 

lower  confidence  Unit  for  p.  Sudakov  (1973)  showed  that  for 

k 

n  <n,,<...<n.  and  g(x)  =  II  (n.-x.), 

1-2-  -  k  i«l 

u(n1,y1,a)  <  b  <  u^.Iy^.a)  , 

k 

where  y1  *  n.^,  qo  “  1  -  n  t(ni~x0 ±) /n^  . 

u(n  ,y  ,a)  is  called  the  Lindstrom-Madden  method  for  determining 
lover  confidence  limits  for  the  reliability  of  series  systems 
(see  Lloyd  and  Lipow  (1962)). 

Lipow  and  Riley  (1959)  used  a  different  ordering  function; 
nevertheless  they  noted  that  for  “small"  n^,  their  tabulated 
values  provided  good  agreement  with  the  results  using  the 
Lindstrom-Madden  method.  For  large  values  of  n^  the  tabulated 
values  that  they  provided  are  based  on  the  Lindstrom-Madden 
method.  Here  we  provide  a  further  justification  for  the 
Lindstrom-Madden  method  By  establishing  that  it  provides  conserva¬ 
tive  lower  confidence  limits  (i.e.  is  a  lower  bound  to  b  defined 
in  (2.9))  using  the  ordering  function  g(x)  employed  by  Sudakov 
and  we  also  obtain  an  upper  bound  for  b,  thus  determining  the 
possible  error  of  the  Lindstrom-Madden  method. 

Sudakov' s  proof  is  unnecessarily  complicated  and  contains 
some  incorrect  assertions,  which  nevertheless  do  not  affect  the 
validity  of  the  conclusion.  In  the  Appendix  we  provide  a  simpler 
proof  of  some  auxiliary  results  needed  for  the  generalization  of 
Sudakov's  theorem  given  below. 
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Theorem  4.1.  Let  g(x)  satisfy  the  hypothesis  of  Theorem  3.1. 
Then , 

*. 


b  <  min  u(n  , [y  .  ] ,o)  , 
l<i<k  1  1 


(4.1) 


where  b  is  given  by  (2.11)  and  y^  «*  y  (j ^ , j ^ > •  •  • , j ^  ^)  is  evaluated 
at  J  ^“0 »  Jt“l ,  2,  .  »  . ,  i*-l.  Note  that  y^^  *  y*.  If  we  also  have 


t  hen 


>  -^1  ,  j-1,2, .  . .  ,k-l  , 

nJ_1j  ~  “j  +  l 


u  ^ni *  yi * 03  - b 


(4.2) 


(4.3) 


Proof :  (4.1)  is  immediate  from  (2.11)  upon  setting  p^  =  l,  j^i  and 

solving  F (xq ; 1 , . . . , 1 , p^ , 1 , . . , , 1)  =  a.  Recall  that  f  n2  1  •••  < 
and 

Iyl3 

F(xo;p)  *  l  b(n1-i1;p1,n1) 


*r° 


Iyk-13 


*•*  l  o0  b(nk-l"ik'-l;pk-l’nk-l)1P.  (nk~tyk3  ’  tyk3+1)  -(4'4) 

*k-l  0  k 

Now,  apply  Lemmas  Al,  A2,  and  A3  to  the  innermost  sum  in  (4.4),  to  get 
tyk-l3 

-  bfnk-l''Ik-l;pk-l*nk-l3lp.  (nk“fyk3  *  tyk3+13  - 
k-1  u  * 

Iyk-13 

J  .  b<nk-l“1k-l8pk-l*Bk-l>1plp<nk"yk*yk+1)  - 
1k-l  0  k 

tyk-l3 

I  b  (n^,!  “l^^i  *Pk-l*nk-l3Ip  ^nk  l""yk  l,yk  l™"^k  l*^*'^'3  ^ 
ik-l"°  k 


1-  -  <nk-l’yk-l,yk-l+1) 


pk-lpk 
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Repeated  applications  of  the  above  establish  that 

F(xo;p)  <  1  k  (n1-y1,y1+l)  .  (4.5) 

n  Pi 

i«l 

(4.3)  follows  immediately  from  (4.5),  completing  the  proof. 

Remarks .  If  (4.3)  holds  and  y^  is  an  integer,  then  b  =  f(n^,y^,a). 

It  has  often  been  suggested  (Lloyd  and  Lipow  (1962), 

Winterbot tom  (1974),  Bolshev  and  Loginov (1966) ,  Mirniy  and 
Solov'yev  (1964))  that  the  confidence  level  should  depend  only  on 
n^,  the  smallest  sample  size.  We  now  provide  a  numerical  illus¬ 
tration  to  show  that  the  bound  in  (4.1)  may  be  improved  by  taking 
all  the  n^'s  into  consideration. 

Let  k=3,  a=.l,  n  ■>  (10,12,30),  x  =  (0,3,0).  Then  for 
3  ° 

g(£)  *  f(n1,Iy1l,a)  -  .541,  f(n2,[y2],a)  -  .525, 

f  ( { y  ^  3  * 01 )  =  .639.  The  use  of  (4.3)  establishes  .500  <  b  <  .525. 

Note  that  if  xQi  **  ni§  for  some  i,  l<i<k,  then  g(x)  =  0  and 
b=0.  It  seems  reasonable  to  use  b=0  as  the  lower  confidence 
limit  whenever  x^^n^  for  any  monotone  ordering  function  satisfy¬ 
ing  the  conditions  of  Section  2. 

k 

We  now  show  that  if  g(x)  =  II  (n  -x  ) ,  then  (4.2)  is  satis- 

i-1  1  1 

fled  and  Theorem  4.1  applies.  This  result  will  extend  a  result 
due  to  Winterbottom  (1974) ,  who  established  this  fact  for 
particular  special  cases.  In  addition,  we  will  also  show  that 
(4.2)  holds  for  a  number  of  other  ordering  functions  used  in  the 
literature. 
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k 

Theorem  4.2.  Let  g(x)  *  II  (n,-x,+a  ),  where  a.>0  and 

i*l  1  1  1 

>  Ct£+^ni»  1*1,2 . k-1.  Then  (4.2)  is  satisfied. 

Proof .  If 


(nj-y.+ot  )  n  (n.+<x.)  =  c 
j-i+1  J  J 


and 


(nj-k^Hn  +o  )  n  (n  +t.  )  .  c 

j  =  i+2  J 


then  we  have 


^ni“yi+ai)  (ni+1+ct1+1) 


(ni-ki+ai) (ni+1-yi+1+ai+1> 


establishing 

yi~ki  =  li± 1  ni-H(lti+aj-ki) 

ni~ki  ni+l  ^ni+l+ai+l^  ^ni~ki* 

Thus  (4.2)  holds  if 


ni+l*ni+Cti~ki) 
^ni+l+ai+l^ni~ki) 


>  1 


this  last  inequality  will  be  true  whenever  ni+iai  >  ai+lni" 
particular,  this  is  valid  when  a^**0,  i  *  l,...,k  which  is 
Sudakov's  ordering  function. 


Theorem  4.3. 

Proof.  If  1 


If  g(x)  ■  1  -  f  x./n  ,  then  (4.2)  is  satisfied 
i-1  1 


-  V"l 


1 - -  - 


i+1 


‘i+3 


,  then 


or 
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This  type  of  ordering  function  has  been  employed  by  Pavlov 
(1973) ,  for  example. 

k  2  i 

Theorem  4.4.  Let  g(x)  *  T  a^x^+z  (a*x.)  ,  where  z  satisfies 
-  °  2^  i  i  a  i  i  a 

l-$(z  )  ■  a  and  $(x)  is  the  standard  normal  distribution  function, 
°  k  . 

>•*•>  *|^»  and  “  (n^  £  1/n^)  .  Then  g(x)  satisfies 

(4.2)  if  and  only  if 


(Wl)yj  -  (aj'aj+l)za+ajkj"“jkj(zaaj+2c'aj(yj+kj))  *  (4-6) 


j-l 


2, 


Proof:  If  g(x_)  **  c  +  £  a^k,,  then  defining  J  ■  c,, 

°  i-1  1  1  i-1  1  A 


Vj  +  WjV*  ■  c 


(4.7) 


and 


”jk3  +  a3+ly3+l  +  2a<cl+alkj+aj+lyj+l)>  ‘  c  •  <4-8) 

Equating  the  left  hand  sides  of  (4.7)  and  (4.8),  we  obtain  (4.6). 
If  k**2 ,  (4.6)  holds  for  all  cases  of  interest. 

If  (4.6)  holds,  then  setting 


,rf(x,l-a)  , 

1-a  -  (r(x))  AJ  tx  Ae  Cdt  , 

0 

a  straightforward  limiting  argument  shows  that 


max  a±f Uyi3+l,l-a)  <  b  <  aLf (yx+l ,  1-a)  .  (4.9) 

This  ordering  function  has  been  used  by  Johns  (1976)  and  b  in 
(4.7)  ir.  the  value  tabulated  by  Johns  for  k*2.  The  validity  of 
the  lower  bound  does  not  depend  on  (4.6).  In  Table  1  below,  the 
lower  and  upper  bounds  given  in  (4.9)  are  tabulated  along  with 
the  values  given  by  Johns  for  a».l.  These  refer  to  upper  confi¬ 
dence  limits  for  the  Poisson  parameter  combinations  ®i^i+a2^2’ 
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Note  in  particular  that  three  of  the  values  tabulated  by 
Johns  (Indicated  by  asterisks)  violate  (4.9).  Specifically  consider 
5.24,  in  which  esse  [y^  -  5,  since  g*(2,5)  -  4.78,  g*(5,0)  *  4.7 2 
and  g  (6,0)  ■  5.48.  Using  the  Poisson  approximation  we  obtain  the 
value  9.275  for  the  upper  confidence  limit  to  X  for  a*.l  and  thus 
al^l+a2^2  "  **.56.  Consequently  the  sup  must  exceed  5.56.  An 
alternative  approach  to  the  one  suggested  by  Johns  for  k  >  3  is  to 
simply  use  a^f (y ^+1 , 1-a)  for  b. 


Table  1 

Comparison  of  Upper  and  Lower  Bounds 
With  Values  Tabulated  by  Johns  for  ct».l 


al 

X1 

x2 

Lower 

Bound 

Upper 

Bound 

Johns' 

Tabled  Value 

.9 

7 

2 

4.79 

5.50 

5.17 

.9 

3 

0 

2.07 

2.27 

2.16 

.75 

6 

3 

6.00 

6.65 

6.23 

.75 

12 

3 

7.90 

8.29 

7.91 

.67 

3 

3 

5.36 

5.61 

5.33* 

.67 

15 

2 

8.71 

9.24 

8.81 

.60 

5 

2 

5.56 

5.62 

5.24* 

.60 

7 

6 

9.24 

9.53 

9.18* 

5. 

Numerical  Examples 

and 

Concluding 

Remarks 

Examples  1  and 

2  Illustrate 

the 

method  we 

have 

described  in 

this 

paper . 

k 

Example  1: 

Let  H(x) 

-  n 

(ni~xi 

)  ,  ct 

«  .05,  k  - 

5, 

n  * 

i-l 

(20,30,40,25 

,60),  x 

-  (2,6 

,10,8, 

15). 

Then  the 

952 

upper  confi- 

dence 

limit 

for  the 

failure  probability  is  contained  in  (.86,. 88) 

3- 


Example  2:  Let  H(x)  «  Il(n^-x^),  a  »  .05,  k  •  2 ,  n  ■  (10,10), 
x  “  (3,2).  Then  the  95%  upper  confidence  limit  for  the  failure 
probability  is  contained  in  (.70,  .73).  The  value  given  in  Lipow 
and  Riley  (1959)  is  .70. 

Remarks ■  In  this  paper  we  have  showed  that  the  Lindstrom-Madden 
technique  is  conservative  for  ordering  functions  satisfying  (4.2). 
Further,  if  y  is  an  integer,  then  the  Lindstrom-Madden  method  is 
exact.  We  have  also  relaxed  the  conditions  needed  in  Winterbottom 
(1974)  and  provided  an  alternative  to  the  method  of  Johns  (1976). 


Append ix 


The  auxiliary  results  employed  in  the  proof  of  Theorem  4.1 
are  provided  here. 

Lemma  Al:  l^(n-x,x+l),  0<y<l,  *s  a  decreasing  function  of  n  and 
an  increasing  function  of  x.  Iy(np,nq+1),  p+q  »  1,  0<p<l,  is  an 
increasing  function  of  q. 

Proof ;  The  proof  is  immediate  from  the  observation  that  the 
beta  distribution  with  parameters  a  and  8  has  monotone  likelihood 
ratio  in  a  and  -6  and  that  if  a  probability  distribution  has 
monotone  likelihood  ratio  in  8,  Fg(x)  is  a  decreasing  function 
of  8  (Lehmann  (1959),  p.  68  and  p.  74). 

^i"^i  ^i+1 

Lemma  A2:  If  - — —  >  -  and  n.  <  n  ,  ,  _  .  then 

n^~ki  -  ni+1  i  -  i+1 

Iy(ni“yi’yi"ki+1)  -  Iy^ni+l~yi+l*yi+l+1)  *  (A.l) 

Proof :  Rewriting  the  left  and  right  hand  sides  of  (A.l)  as 

S  <ni'ki)(1  -  i 

”l+l(1  -  •  <*•*> 


Lemma  Al  applies  and  the  conclusion  follows. 
Lemma  A3:  Let  y^  -  y,  (Ky^l,  i-1,2.  Then 


Jv  v  (n-x»*+D  >  1  *>(n-k;y,  ,n,)I.  (n-x,x-k+l)  .  (A. 3) 

2 


Proof : 


.n.  n-k..  .k  Kn-k+l)  fy2  ^n-x-1,,  ^x-kJk 

I  <k>'i  r(n-»)r(»-k+T)  '  <l  °  dt 

k**0 


k  n-k 

k“0 

y,y0  1*]  (l-y.)~t 


r(n-t-l)  (1~yi)  yi  fyLy2  ft]  n-x-1  f  yl~t\  *~k  .dt 

r(n-x)  .£0  k!  T  (x-k+1)  J  XyJ  V  y1  /  yx 


r(n+l)  fyly2  tx~yl 

J  k«0  kl 


0 

k  n-x-1 


(y^t) 


x-k 


r (n-x) 


ru-k+1) 


dt 


Thus  (A. 3)  will  hold  whenever 


T<n 


£1^  pT*  » 


r  i  /i  \ k  n-x-1 ,  ..x-k 

r<n-m  Piyi  (?]  (1-yi)  1  4t 

r (n-x)  J  k£0  k!  T(x-k+l) 


or 


T(n-H)  fy,y 


r (n-x) r (x+i) 


yly2  tn-x-l(1_t): 


r  (*>  r(,.ki)  >  o  . 

I  k-o  k’r<*-k+i)  U-t  /  Vi-t  ;  J 


k-0 
Ti-* 


y,-t  /  i-y,\ 

Writing  j—-  •  (l  -  j, -_-t  ■  J  and  noting  that  0<t<yjy2<l,  and 

(l-t)>l-y1»  ve  observe  that  (A. 4)  holds  and  the  lemma  is  proved. 


(A. 4) 
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